function dl = ra_to_dl(R, a, B)
dl = a*(2*R-B)/2;
end